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In this paper we provide a consistency relation between the amplitude of the hemispherical bipolar 
asymmetry, A, and the amplitude of the primordial non-Gaussianity in the squeezed limit, /jvi, as 
MM ~ 10 _1 /jvl- We demonstrate that this consistency condition is at work for any model of 
inflation in which the curvature perturbations is sourced by a single light field with the Bunch-Davies 
f—*) . initial condition, irrespective of the number of inflation fields which contribute to the background 

^v^j ' inflationary expansion. As a non-trivial example, we show that observable hemispherical asymmetry 

can be generated in single field non-attractor inflationary models. We also study hemispherical 
asymmetry generated in the models of multiple fields inflation. We show that A is controlled by the 
weighted sum of non-Gaussianity contribution from each field. In particular we show that observable 
hemispherical asymmetry can be generated in models where inhomogeneities are generated from a 
light scalar field modulating the surface of end of inflation. 
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I. INTRODUCTION 



o 

u 

Cosmology is in the era of precision. The release of the Planck data [1[ confirmed the validity of the Standard 
Model of Cosmology, ACDM, with the initial conditions from inflationary paradigm Q seeding the late time Large 
| Scale Structures (LSS) of the Universe, which are nearly Gaussian, adiabatic and scale- invariant @. However, the 
■ Planck collaboration also reported a hemispherical asymmetry in Cosmic Microwave Background temperature (CMB) 
fluctuations [|[ which was also observed in Wilkinson Microwave Anisotropy Probe (WMAP) data The statistical 
significance of this finding may be under debate 0, 0- However, should the upcoming observations confirm this 
finding, then this may be viewed as a challenge for simple models of inflation which predict an isotropic Universe with 
(*T) . all pre-inflationary history washed out during the near exponential expansion of the background (all. 
t-H ' The observed asymmetry may indicate to a pre-inflationary physics. Historically Grishchuk and Zerdovich[Ioj 
first considered the temperature anisotropy from superhorizon perturbations, then Erickcek et al. [TT| showed that 
a super- horizon perturbation may introduce this asymmetry, the idea which also mentioned in Gordon[l2l|. This is 
l/"} \ because the perturbations produced by the inflaton field depend on the background value of the field. Therefore the 
i long wavelength mode (i.e. a super-horizon mode) changes the perturbations and are potentially capable to introduce 
asymmetry in curvature perturbation power spectrum. The decomposition of modes into the long (super-horizon) and 
short (CMB) scale modes and their add up is very similar to the Peak- Background Splitting of the modes in LSS [Hj], 
where, in analogy, the long wavelength mode (horizon size mode) changes the amplitude of the matter perturbation 
on short scales (structures' scales "i.e. cluster of galaxies"). 

Recently a very interesting proposal was put forward by Schmidt and Hui [l4[, that the modulation from the 
primordial non-Gaussianities (PNG) can produce the asymmetry by mixing the long and short wavelength modes [l5j , 
for a related idea see also |l6[. This is an interesting observation which can yield a bridge between the predictions 
from PNG and observations. However, it is worth to mention that the asymmetry generated from the modulation 
of the large amplitude long wavelength mode will introduce large angular temperature fluctuations, like quadrupolc 
and octupole anisotropics on CMB, so the observed uniformity of CMB map constrain these departures from isot ropy 
17 1 . These constrains mainly come from the quadrupole, C2, and the octupole, C3, moments of CMB [ToL ll7T - fl9l |. 



This is because the primordial curvature perturbations are translated into the perturbations in gravitational potential 
($) which causes the temperature fluctuations by Sachs- Wolfe (SW) effect (AT/T ~ 3>/3) [2(|. This gravitational 
potential modulation can not generate a monopole and a dipolar in CMB power spectrum. This is because they are 
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canceled by the angle avera ged CMB power spectrum and the intrinsic dipolar term introduced from the peculiar 
velocity respectively 0, EHTHI- However, the quadrupole and the octupole terms are distinguishable and also are 
well constrained by CMB observations [22j . 

A phenomenological parametcrizations of the dipolar asymmetry is defined via 

V][ 2 {k^) =V l 1 l 2 "\k){\ + A(k)p.x/x n ) (1) 

1 /2 iso 

in which Vn(k,yL) is the curvature perturbations power spectrum, is the isotropic power spectrum, A(k) is 

the amplitude of the bipolar asymmetry, p is its direction and x n is the comoving distance to the surface of last 
scattering. The Planck data indicates A = 0.072 ± 0.022 on large angular scales, £ < 64, with the best fit for the 
anisotropy direction (L b) = (227,-27) [|. 

In Erickcek et al. it is shown that the maximum anisotropy generated in models of single clock inflation 

is A max ~ 0.01 in order not to contradict with the observational constraints from the quadrupole and octupole 
constraints. This value of A is almost an order of magnitude smaller than the central value of anisotropy amplitude 
reported by Planck. One can easily check that in conventional models of single filed inflation A oc (n s — 1) in which n s 
is the curvature perturbation spectral index, which from the Planck observation, is obtained to be very close to unity. 
As a result, in simple single field inflation models, the smallncss of A is directly related to the smallness of n s — 1. If 
one compares this situation with the Maldacena's consistency relation [23| for the amplitude of PNG in the squeezed 
limit, /atl, in which f^L ~ n s — 1, the similarity rings the bell. To make this similarity more pronounced, it is observed 
that in curvaton models with large enough value of Jnl an observable value of A can be generated [Tl|, [25[ . These 
results tempt to indicate that in general the amplitude of A is controlled by Jnl- This was directly demonstrated by 
Lyth in curvaton scenario in which he found A oc Jnl- We prove this conclusion in a model- independent way. We 
show that for all models of inflation in which the curvature perturbations have a single source with the Bunch-Davies 
initial condition this consistency condition does hold. 

There are additional constraints from LSS on the asymmetry from power spectrum modulation. The SDSS sample 
of quasars sets a very strict constrain on A in the scales of k ~ \Mpc~ 1 [27} ■ Although the hemispherical asymmetry 
shows up in very large scales, k < 0.035Afpc _1 , but this constraint implies that the power modulation must vanish 
very quickly in going from large to small scales. As a result, the two cosmological observations: a) CMB anisotropy 
constraints from quadrupole and octupole and b) the constrains from the LSS (i.e. quasars) on the primordial 
anisotropic power spectrum in small scales, restrict the space of the inflationary models based on super-horizon 
modulations of V-R,(k, x). On the other hand, this may open up a new venue to investigate other cosmological models 
as the source of asymmetry, like considering the modulations of the spectral index, the modulation of the tensor 
perturbation amplitude, the modulation of the optical depth or spontaneous symmetry breaking [8l. [28j. 

Note added: While this work was in its final stage, the paper [2!| appeared which has some overlaps with our 
multiple fields analysis in Section [TTT1 



II. SINGLE SOURCE INFLATION AND THE CONSISTENCY RELATION 



In this section we prove our consistency condition for models of inflation in which the curvature perturbation is 
generated by a single source. Of course there may be more than one inflaton fields at the background which can 
contribute to background expansion to solve the flatness and the horizon problem. But we assume that only one 
field is responsible for perturbations in order to generate the observed CMB fluctuations and the structure formation. 
With this definition, the standard curvaton model [3(1 HH is also considered as a single source scenario in which the 
curvature perturbations originate entirely from the curvaton fluctuations. On the other hand, if both the inflaton 
field and the curvaton field contribute to perturbations, then it is a multiple source scenario which is studied in next 
section. 

Having this said, one may wonder what is so crucial about the assumption of single source inflation. To see this, let 
us consider the curvature perturbation on comoving surface 1Z. With an appropriate choice of gauge transformation, 
one can always find a coordinate system (comoving gauge) in which the comoving velocity potential Su vanishes. 
Furthermore, in models in which perturbations are sourced by the single field 0, one has Su oc S(f>. Now, in the 
comoving gauge in which Su = 8<p = 0, the re-normalized scale factor a(i,x) on each Hubble patch (in the sense of 
separate Universe approach) is given by 

5(t,x) = a(i)e TC(t < x) (2) 
in which a(t) is the background scale factor. As a result, the three-dimensional spatial metric on the surface of 
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constant time is given by 

ds 2 {3) = a(t) 2 e 2TC( *< x Wd Xl . (3) 

In this view, the effect of a large scale curvature perturbations will be a modulation of comoving distance x — > e TC ^ ,x ^x 
which will play essential roles in our analysis below. 

Our goal is to see the effects of very large super-horizon scale perturbations on smaller CMB scales perturbations. 
In this view, the size of our observed Universe is at the order of Hq 1 in which Hq is the current Hubble expansion 
rate while the long mode (the modulating mode) has the wavelength Xl ^> H^ 1 . As usual, in order to perform the 
calculations we go to the Fourier space with the volume V . The volume of Fourier space V should be larger than the 
whole universe. On the other hand, we also assume that the scale of long wavelength mode is much larger than the 
scale of Fourier space. As a result we have the following hierarchy 

H - 3 « V « (4) 

in which = 1/Al is the long mode comoving wave number. As a result, the Fourier transform of long wavelength 
mode satisfies the following relation 

ll L ( X )~K kL e^* (5) 

It is worth to mention that we have assumed one specific mode fcz, with a known wavelength in real space. Consequently 
there is a delta function term in Fourier space integral 8{k — Kl), resulting in Eq. ([5]). Now our aim is to parameterize 
the effects of fcz, on CMB scale modes k such that 

TZ = TZ k (1 + A(k)p.x/x n ) , (6) 

which is equivalent to Eq. ([1}. Note that curvature perturbation on logarithmic momentum scale IZk is defined via 

(K k K k ,) = (2TT) 3 6 3 (k + k')P K (k) , V K = ^P n (k). (7) 

Starting with the parameterizations Eq. ([1]) the fractional change due to long wavelength modulation on curvature 
perturbation power spectrum is 

which will be used later. 

Now we can study how the local non-Gaussianity modifies the power spectrum of CMB-scale perturbations, leading 
to anisotropy. To start, let us assume that the curvature perturbations does not evolve on super-horizon scales, i.e. 
IZk = for k <C aH . This is true for all single field models in which the system has reached the attractor regime in 
which the decaying mode can be neglected. However, this assumption does not hold in models of non-attractor single 
field inflation H, [33| in which the would-be decaying mode actually dominates over the would-be growing mode, the 



constant mode. As a result, during the non-attractor regime TZk ^ even for single field inflation. We will study the 
general non-attractor case shortly. Consider the three-point correlation function of the long wavelength mode with 
two small (CMB scales) modes as follows 

n^M^n^)) * ( k^) ( k^k^J * (^(k L) (^(k L )^y(^(k 1 )7e(k 2 ))| TC(kL) ^) ) ( 9) 

This equation is the key for our analysis below. Intuitively this equation means that the effects of a large scale 
perturbation would be just a rescaling of the background for small scale perturbations. This is because the small scale 
perturbations can not probe the spatial variation associated with the long wavelength mode. This logic was employed 
in [33 - [36| to calculate the amplitude of the local non-Gaussianities and to check the Maldacena's consistency condition 
that Jnl^ (n s — 1)- Note that to write Eq. ^ we have assumed that there is only one degree of freedom encoded in 
TZ, so there is no other source of perturbation. Secondly, it is assumed that there is no intrinsic non-Gaussianity deep 
inside the horizon which means we have started from a Bunch-Davies vacuum. In other words, the non-Gaussianities 
generated in the model are either at the time of horizon crossing or on the super-horizon scales. These two assumptions 
(the Bunch-Davies (BD) initial condition and the assumption of single source perturbations) are the only assumptions 
in our analysis in this section. However, note that we do not impose the slow-roll conditions. As long as a violation 
of slow-roll condition does not violate the BD initial conditions for modes inside the horizon, then our consistency 
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condition docs hold. For example, this is the case when the violation of slow-roll conditions happens after the modes 
of interest has left the horizon. 

Note that Eq. ([9]) may not be a practically useful method to calculate /nl- We instead employ an inverse 
engineering method in which we look for the amplitude of asymmetry, A, for a given value of Jnl- Therefore, to 
calculate the bispectrum one has to employ other methods, such as the in-in formalism or the 6N formalism. 

Defining the bispectrum i?7j(ki, k 2 , k^) in the Fourier space via 

(1l ki n k2 n kL ) = (2 7 r) 3 <5 3 (k 1 +k 2 +k L )B K (k 1 ,k 2 ,k L ) (10) 

Eq. © yields 

B R (k 1> k 2> k i )=J^ fcl ^. (11) 
On the other hand, the amplitude of the local non-Gaussianity f^L in the squeezed limit k^ <C k\ ~ fc 2 is given by 

fNL = i im ^B n {^MM) ( ) 



Comparing Eqs. (|TT|) and (fl2|) yields 

1 dP nk _ 12 
Pn k OTIl 5 



fNL (13) 



in which it is understood that k = k\ = fc 2 corresponds to the CMB-scale perturbations. 

Since we are looking for the directional modification of power spectrum due to large scale perturbations we can 
write the above equation as 

VP TCfc = T^VKi • (14) 



Now comparing this equation with Eq. © yields the following formula for A 

A = lf NL x n \Vn L \ (15) 
o 

On the other hand, using Eq. (J5)) for the decomposition of the long mode in the Fourier space, we have 

\viz L \ ~ k L \n L \ = k L vUl = k LEPn 2 » (is) 

in which we have used TZl = and following the Lyth convention p5| . we have introduced the parameter E via 

Pn L = E 2 V n . (17) 

Here it is understood that 1Z is the CMB-scale curvature perturbation with the power spectrum V-ji- Now plugging 
the value of |V7£l| obtained in Eq. (|T5)) into Eq. (|15p yield our desired formula 

A(k) = ^f NL k L x n ET^ (18) 

This is our consistency condition, relating the amplitude of bipolar hemispherical asymmetry to the amplitude of 
non-Gaussianity in the squeezed limit. Note that this formula applies for all single source models of inflation with the 
BD initial conditions which have reached the attractor limit, i.e. with the assumption 72. = on the super-horizon 
scales. Note that Eq. (fT5)) was first obtained by Lyth [25| for the curvaton scenario. 

Having obtained our consistency relation for the attractor models, now we proceed with the analysis for the general 
single source inflationary models including the non-attractor models as studied in [321 . |33| . The analysis go parallel 
to the above analysis, except that in our starting point, Eq. ©, we have to also consider the modulation of the 
two-point functions of small scales by 1Z. Assuming that both 1Z and TZl are independent variables, and following the 
same procedures as before, yields 

(^(k L )^(k 1 )^(k 2 )) ~ (ft(k L ) ^(k L )^^(^(k 1 )7e(k 2 )) + 7e(k L )^^(7e(k 1 )^(k 2 ))) ) (19) 
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As before, using the definition of the bispectrum from Eq. (|10j) . this equation can be further manipulated to give 
12 f p p 1^ (v dP K ^nb 9Pn\\ p 5P K 1 dP n , 



Now we need 



r) Pn f) Pt> 

vp —wt^ + w; VKl (21) 



1 /2 

Noting that V72l = khT^L and TZl = T > ii L we conclude 



Comparing this with ([2U]) yields 



^ - ^fNLkLVUl = l 4fNLVn . (23) 



P n ~ 5 JJV " " 5 



Interestingly, this is exactly the same as (|14p and following the same steps as before, we obtain our consistency relation 
Eq. (TTgj) . Note that in above calculations, in going from Eq. (f2"Tj) to Eq. (f2"3"]l . we have used the fact that 72 z, and 72 z, 
are in phase, i.e. 72l(x) = 1Zk L {t)e~ L '* and 72z,(x) = TZk L (t)e~ lkLX which follows from our long wavelength mode 
decomposition given in Eq. (|5|). 

In summary, we have the consistency relation given by Eq. (|18|) valid for all models of inflation in which the 
perturbations are produced by a single source with the BD initial condition. This result is independent of the slow- 
roll assumptions and also is valid for non-attractor models alike. Furthermore, this result is independent of the 
number of background fields which can contribute to the background inflationary expansion but do not contribute to 
the curvature perturbations. This include the standard curvaton scenario in which the inflaton docs not contribute 
to the curvature perturbations. 

To obtain useful information from the abstract formula Eq. (|18[) we need to have some information on the mag- 

1/2 

nitude of the combination k^XnEV^ . The Observational constraints from the quadrupole and octupole moments 
of anisotropic power spectrum of CMB can be used to constrain this combination. To do this, note that we are 

1 /2 

considering a specific mode of long wavelength modulation (i.e. Tic = PP-jz sin(fcx,.x + coj) and its gradient in the 

CMB is V72l — \(kx n )EV^ 2 cos(w)| with uj being an arbitrary phase. In Erickcek et al. pj} it was shown that the 
quadrupole and octupole constraints are as below: 

{kxnfE-p^ 2 sinw < (5/3) x 5.8Q 2 - (5/3) x 5.8 x (6 x 1(T 6 ) (24) 
{kx n fET]l 2 cos u < (5/3) x 32Q 3 - (5/3) x 32 x (9 x 1(T 6 ) 

where the additional factor 5 /3 comes from exchanging the gravitational potential $ to curvature perturbation 72.. We 
note that the phase appearing in quadrupole and octupole are the same as in the long wavelength mode. The octupole 
can constrain the change of the long wavelength mode independent of the phase, while the quadrupole vanishes for 
uj = 0. Consequently, the constraints on the power of long wavelength comes from octupole term. On the other, hand 

1/2 

the amplitude of long mode perturbation inside the big box, i.e. the super-universe, is 72l = EV^ . In order for this 
to be consistent as a perturbations, we require 

72| ~ E 2 V n (ki) < 1. (25) 
Combining the bounds from Eqs. (j2"4")l and (f2"5"j) yields 

\{k L x n ) EV][ 2 < 32Q 3 1/3 - KT 1 (26) 

Plugging this upper bound into Eq. (|18j) , we obtain the following upper bound consistency condition for the amplitude 
of dipolar asymmetry 

\A\<lO^\f N L\. (27) 
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To explain the observed anisotropy we require |^4(A;)| = 0.07 ± 0.02. Now we employ our result to some interesting 
examples. 

First consider the standard slow- roll single field inflation. From Maldacena's analysis [H| (see also [3]) it is known 
that f NL ~ (n g - 1). Using this in Eq. (J27J) gives A < 10 _1 (n s - 1) - 10~ 3 , as observed in [H[25|]( see also [26| in 
different context) . This is too small to explain the observed dipolar asymmetry. Now consider the standard curvaton 
scenario, in which all perturbations are sourced by the curvaton field. In this model, Jnl is uncorrelated to n s — 1 
so one can easily saturate the observational bound from Planck Jnl = 2.7 ± 5.8 (68 % CL) [3{|. As a result one can 
easily saturate the observational bound |A(/c)| = 0.07 ± 0.02 on the CMB scale by taking |/jvl| ;$ 5 (25|. However, 
/at l is nearly scale-invariant for the standard curvaton scenario so this model may not work easily in order to reduce 
the amplitude of asymmetry on small scales as required from the quasar constraints. 

As a non-trivial example, we consider the non-attractor single field models studied in [12, Hj|. In these models, 
inflation has two stages: the non-attractor phase at the early stage of inflation followed by the attractor regime 
towards the end of inflation. During the non-attractor phase the curvature perturbation is not frozen so 1Z is not zero 
on super-horizon scales. As observed in [H, HH a large value of local non-Gaussianity can be generated given by 

_ 5(1 + eg) 

■> NL ~ £3 ' ( 28 ) 

in which c s is the sound speed of perturbations during the non-attractor phase. Choosing small enough c s one can 
easily saturate the observational bound on /jvx. For example with /jvl ^ 2 one can easily find large observable bipolar 
asymmetry. One interesting aspect of this model is that for the modes which leave the horizon during the attractor 
phase the curvature perturbation is frozen so 72. ~ once these modes leave the horizon. This means that /jvl has 
a built-in scale-dependence in non-attractor scenarios. For modes which leave the horizon during the non-attractor 
phase fi^L is given by Eq. (|28[) while for the modes leaving the horizon during the attractor phase /jvx — n s — 1 ~ 0. 
This helps to satisfy both the CMB constraints and the quasar constraints on A(k). For this picture to work, we tune 
the parameters such that the observed CMB scales leave the horizon during the first 5-10 numbers of e-folds when the 
system is in non-attractor phase. Scales smaller than Mpc arc supposed to leave the horizon once the background 
reaches the attractor regime during last 50 or so e-folds. By construction /jvl an d consequently A{k) are reduced 
sharply for these scales. It will be an interesting exercise to fit the predictions of this scenario with the Planck data. 



A. SN formalism for single source models 



As a further illustration, here we employ the SN formalism [4CJ, |4l| to obtain Eq. (|T8l) for a generic single source 
inflationary model. Suppose the perturbations are generated by the single source (in the sense defined above) Sip. We 
have to Extend the usual 8N formalism to more general case in which the background number of e-fold N is defined 
in the phase space [12, HH so N = N(<p, <j)). Note that the dependence of TV to p> is crucial in non-attractor models as 
emphasized in [3^ , otherwise one obtains a wrong result for 1Z. Expanding N to second order in perturbations of 8(p 
and S(p, the comoving curvature perturbation to second order is given by 

K = + N^54> + \n^8<P 2 + \n m 5^ + N^8(f)Sj> . (29) 

In this view both 8<p and 8(p are treated independent random variables. However, in separate Universe approach 8<p 
and 8(f> are classical fields measuring the deviation from a given background solution. In the gradient expansion limit 
where both 8<f> and 8<f> becomes classical, one expect that 8(p carries the same initial quantum (statistical) information 
as 8(f). As a result, to leading order in gradient expansion 8<p and 8<p are related to each other via 8<f> = a{t)8<p in 
which a(t) is a function depending only on background quantities. For example in models of slow-roll inflation in 
which 8(f) freezes out on super-horizon scales a = 0. On the other hand, for the non-attractor model studied in [IH 
8(f) = 8(f>oe aN with a a constant. With this assumption in mind, from Eq. (|29|) one obtains 

6 , + 2aN M + a>N M 

tJnl = -2 • (30) 



Note that to obtain this formula, we do not to know the specific form of a(t). All we need is that a(t) does not carry 
independent quantum (statistical) information. 
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On the other hand, expanding the coefficients of linear parts of TZ in Eq. (|29|) as a function of background field 
modulated by the long mode S(j)L, i.e. N^((f> + 8<f) L ) = N^{(f>) + Sc^lN^, one obtains 



K ~ TV 



1 



1 + pNL (Nj + aN ^ 5<j>! 



= K lso 

Following the decomposition of long wavelength mode in Eq. ([5]) we have 

(50l(x) = V^, 

Plugging this into Eq. (|3~Tj) yields 

)l/2 iki,.x 



Vn = V 



■R 



As a result, the fractional change in the gradient is 

V7> K 12 



Vn 



1 /2 

As before, we have to relate the power spectrum of the long mode, VA , to the power spectrum of the short mode, 



(31) 



(32) 



(33) 



(34) 



-r.1/2 . -r>l/2 

comparing with the definition given in Eq. ([8]) yields 



1 /2 _____ 

EVA . Plugging this relation in Eq. (|3"4"|) . using the identity V t 



1/2 



^ 2 /(N^ + aNtf and 



A{k) = -f NL k L x n EV%*. 



(35) 



Interestingly, this is exactly the same as Eq. ([18 



III. ASYMMETRY FROM MULTIPLE FIELDS MODELS 



In this section we study the asymmetry generated in models of multiple fields inflation in which more than one 
field contributes in curvature perturbations. As mentioned before, the example of curvaton was studied in the 
literature and it is argued that the asymmetry generated from curvaton scenarios with a mixture of inflaton and 
curvaton contributions to the curvature perturbations can generate observable bipolar asymmetry fTA [25l . _4__J . Here 
we consider the general case of multiple fields scenarios in which the quantum fluctuations of more than one light 
scalar fields source the power spectrum and the bispectrum. The goal is to see if the super-horizon modulations 
of these light scalar fields can generate hemispherical asymmetry and to find the relation between the amplitude of 
bispectrum and the amplitude of the asymmetry. For a related recent work see [4__J. 

To be specific, we consider the model containing two light scalar fields (f> and a which contribute into both background 
and the perturbations dynamics. The generalization to more than two fields is straightforward. Employing the 
standard 5N formalism for super-horizon perturbations, the curvature perturbation is given by 

K = N^6<f> + N„5a + \n^54> 2 + ^N^Sa 2 + N^S&a (36) 

To simplify the analysis, we assume that N does not depend to (f> and a so the system has reached the attractor 
regime. The power spectrum is given by 

Vn = N%V S ^ + N 2 a V Sa , (37) 

in which Vsrp and Vga represents the power spectrum of light fields quantum fluctuations Vs<p = Vga = (H/2tt) 2 . 
However, one may consider the gen eral case in which the light scalar fields do not have equal power spectrum. 

Following the notation of [iil. |45| (see also [f___|), the non-Gaussianity can be written as the weighted sum of each 
field's contribution 

fNL = w\j% L + w\r N L + ^W<f w <rfNL ( 38 ) 
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where, we have defined 



6 



G 



G 



5 N% ' 5 JfTL ~ ^ 

and w a and are the fractional contribution of fields a and 



N r 



to the power spectrum, respectively. 



N 2 



N 2 + Nl 



1 



W<j, . 



(39) 



(40) 



where the last equality in holds if we assume = Vs a = {H/2ir) 2 . 

Now we can obtain the anisotropy of the model. Assuming a modulation of the background by two large scale 
modes 8<jl and <5</>l, one has 



K = N^8(j> + N, a Sa + N^S^S^l + N^ a SaSa L + N.^S^Sol + 54> L Sa) 
For large scale perturbations we have 

HL{x)=V] / (p 2 L e'^ , Sa L (x) 
As before, we can relate the power spectrum of the large scale to the power spectrum of small scale via 



.pl/2 ;k L .x 
' Scr L c 



— E' 



V 5 a L = E 2 Vsa 



(41) 



(42) 



(43) 



where we have assumed that the enhancement factors are the same for both fields. This seems a reasonable assumption 
because the dynamics of generating the light scalar fields quantum fluctuations deep inside the horizon from an initial 
Bunch-Davies vacuum and the subsequent horizon crossing effects are the same for both light fields. However, one 
may consider the general case in which the enhancement factor E amy be different for the two scalar fields. Therefore, 
we have 



K = 



l/2\ _ik L 



E ( N^VH 2 



4><* r &<t> 



l/2\ -ik L 



5a 



Using this form of 1Z, one can check that the directional dependence of Viz is given by 



= x 1ev. 



1/2 
■R 



N., 



(44) 



(45) 



where sgn is the sign funtion, +1 for postive argument and —1 for negative. Comparing this with the anisotropic 
ansatz for power spectrum asymmetry Eq. ([8]) yields 



.4 



f^iwaw 1 ^ sgn(^) + wl /2 W< p) + f^w 3 / 2 sgn(^) + f^ L w\ 



N, 



3/2 



N r 



< W 



G 



(46) 



in which we have used the octupole constraint Eq. (|26l) to eliminate the factor — x n kLEV^l 2 . This is a very interesting 

result which is also intuitively reasonable. The amplitude of asymmetry depends on how much each field contributes 
to the curvature perturbation power spectrum, parameterized by and w a , and also on how strong the coupling 
between the small modes to the large modes for each field's perturbations is, as parametrized by the corresponding 
/jvl- Interestingly, we also see the effects of coupling of small modes to the large modes with two different sources, as 
indicated by the term f^ L in Eq. (|46|) . Note that our result Eq. (|46|) is different than the result obtained by Erickcek 
et al. in the particular case of curvaton in which they found A oc fjj 2 . 

As discussed in [29[ one can chose the signs of ffj L , f^L^ ftfL an< ^ the weights of each field in power spectrum such 
that the total Jml in Eq. (|38| can be small enough to satisfy the Planck observations Jnl = 2.7 ± 5.8 (68 % CL) 
[39j] . Then one still has enough free parameters in Eq. (|46|) to make A large enough to saturate the observational 
bounds deduced from Planck. On top of this, one should also check the observational constraints on the amplitude 
of tri-spectrum tjv l which in this model is calculated to be 44 1 



25 
3G 



2 fcr<fi f<t> 



(47) 
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The upper bound inferred from Planck data is [39| tatl < 2800 (95 % CL). It would be an interesting exercise to 
search for parameter values in which large observable value of A can be obtained from Eq. (|46p while the Planck 
constraints on /jvl and tjvzo given respectively by Eq. (|38|) and Eq. (|47|) . arc satisfied. 



As an interesting example, consider the case in which the inflaton field 4> is Gaussian so fj^ L ~ f^ L ~ 0. From Eq. 
38l) one obtains f^ L ~ Jnl/w 2 so Eq. (|46|) yields 



^<tt^=- (48) 
10\/uw 



This is a very interesting result. If one reduces the contribution of the non-Gaussian field a in the power spectrum, 
i.e. w a <C 1, with /jvl at the order of unity one can get large upper bound on the amplitude of asymmetry. 
It may be more transparent to eliminate w a in terms of t^l in Eq. (|48p . Using Eq. (|47p . one obatins 



A < . (49) 

~ 12 v 1 

It is very interesting that the upper bound on A is independently controlled by tnl ■ With the upper bound tjvl < 2800 
(95 % CL) from Planck data |39| one obtains A < 4.4. 



IV. ASYMMETRY FROM INHOMOGENEITIES GENERATED AT THE END OF INFLATION 

As a concrete example, in this section we study bipolar asymmetry generated from inhomogeneities generated at the 
surface of end of inflation. This model was studied in [47H49| as a mechanism to provide large local non-Gaussianity. 
Here we follow the model studied in [i^ and [50} . 

The model is based on inhomogeneities generated from the modulation of the surface of end of inflation in the 
background of a hybrid inflation [5lT l52l| scenario. The model contains three fields, the inflaton field <j>, the waterfall 
field x an d the light isocurvature field a. The potential is 



V 



X 



X 



A 



a 2 -v 2 

\<t>\ 2 + \x 2 ° 2 



(50) 



In this picture, inflation proceeds as in usual hybrid inflation models. The waterfall field \ is very heavy so it quickly 
rolls down to its local minimum \ = 0) so during the whole period of inflation we set x = 0. Inflation ends when the 
inflaton field reaches the surface of end of inflation defined by 



72 „ 2 

7 a 



in which 4>„ indicates the value of 



at the surface of end of inflation and the critical value 

A - M 
9 



is defined by 



(51) 



(52) 



Once 4> hits the surface of end of inflation, the waterfall fields becomes tachyonic terminating inflation very quickly. As 
demonstrated in recent literature (53l - [60j the waterfall quantum fluctuations do not contribute to large scale curvature 
perturbations. 

The idea in [47M49I ] is to generated inhomogeneities from 5a quantum fluctuations at the surface of end of inflation 
which also translate into 8(f> perturbations in Eq. ([ST]) . With these descriptions, the final result for the curvature 
perturbation at the end of inflation, £ e , is the following (see [5(| for details). 



Ce 



350* 3 7 2 er 



ag 



2 J.2 



6(T 



I 



1 — J- a 



T 



2a g^l 
5a F{1 -i 



1 



2(1 -FY 



ya 2 



2^2 
2' 



5a 1 



(53) 



in which we have defined 



a = 



H 2 



Yl\m 2 Ml 



(54) 
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and 



T = 



l 2 °l 



(55) 



Note that in Eq. ([53)) we have neglected the second order term from 5(f> 2 since the non-Gaussianities generated by 
inflaton is very small. Also note that a measures the mass of the infiaton compared to H and as usual we assume 
that a <C 1 in order to get long enough period of inflation. Furthermore, the parameter T defined above measures 
the contribution of a in energy density at the end of inflation, see [5(| for details, so < T < 1. 

Assuming that the two fields S(f> and 8a are canonically normalized with the same amplitude, the contribution of 
the modulator to the power spectrum defined by Eq. (|40[) is here given by 



T 2 cf>l + (1 - F) 2 
and the resulting non-Gaussianity, using Eq. (|38[) is the following 

\f NL = \wlr NL (57) 

with 

V NL (58) 

For fixed w a , and J< 1 wc obtain a large non-Gaussianity. 

Plugging the above value of $nl in Eq. (|46p , the upper bound on the amplitude of asymmetry is 

|A|<10-x w 3/ 2 ^|±Z) (59) 

So we can obtain large asymmetry by requiring J- <C 1. Note that, using (|58[) . the last equation can be rewritten by 

\A\ < . (60) 

1 1 ~ 10 

This is identical to Eq. (|4"5]l which is expected since in this model the inflaton field is Gaussian and the non- 
Guassianity is generated from the light field a by modulating the surface of the end of inflation. Similarly, eliminating 
w„ in Eq. ([60]) in favor of tjvl yields Eq. (|49|) . As described below Eq. (|49p . large upper bound on A is obtained 
with the current Planck constraints on tjvl- 



To summarize, the observed large scale anomalies in CMB may indicate to a pre-inflationary physics. In this 
work we have shown that a long wavelength super-horizon mode with a large amplitude can generate anisotropies on 
CMB scales, which can be related to the local type non-Gaussianity. In this work we have found an upper bound 
consistency relation between the anisotropy amplitude seen by Planck and the non-Gaussianity obtained from single 
source inflationary models. We have used the bounds from octupole moment of CMB as a crosscheck to find the 
upper bound consistency relation. We have shown that a large upper bound on A is obtained in single field non- 
attractor models of inflation which is consistent with the observational constraints from CMB and LSS (i.e. quasars) . 
This is because in these models the usual consistency condition /jvl ~ 1 — n s does not hold and large observable 
non-Gaussianity can be generated. 

We also studied the asymmetry generated in models of multiple fields inflation in which more than one light field 
contribute to the curvature perturbation. It is shown that the amplitude of hemispherical asymmetry is controlled by 
the weighted sum of non-Gaussianity generated by each field. In the special case in which only one field contributes to 
non-Gaussianity while the other field is Gaussian, the amplitude of bipolar asymmetry is controlled by tjvl as shown 
in Eq. (|49|) . In particular, if one reduces the contribution of the non-Gaussian field in the curvature perturbation a 
large upper bound on A is obtained as shown in Eq. (|4"5)l . As a concrete example of models of this type we studied 
asymmetry generated from inhomogeneities from the modulation of the surface of end of inflation. We have found 
that one can generate an observable asymmetry with /jvl ~ 0(1) if the contribution of the modulator field to the 
power spectrum, w a , is sufficiently small. 
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